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Abstract 

A two-component fluid representing dark energy is studied. One of 
the components has a polytropic form, while the other has a barotropic 
form. Exact solutions are obtained and the cosmological parameters are 
constrained using supernova type la data. In general, an open universe 
is predicted. A big rip scenario is largely preferred, but the dispersion in 
the parameter space is very high. Hence, even if scenarios without future 
singularities can not be excluded with the allowed range of parameters, a 
phantom cosmology, with an open spatial section, is a general prediction 
of the model. For a wide range of the equation of state parameters there 
is an asymptotic de Sitter phase. 

Pacs numbers: 98.8G.-k, 98.80.Es 

1 Introduction 

Several cosmological observables indicate that the present Universe is in a state 
of accelerated expansion. The first evidence in this sense came in the end of 
the last decade, when two independent observational projects [1], using type 
la supernovae luminosity distance-redshift relation, provided an estimation of 
the deceleration parameter q — —aa/<P. With a catalogue of about 50 type 
la supernova with low and high redshifts, the analysis revealed a negative g, 
indicating an accelerated Universe. Today, more than 300 type la supernovae 
have been identified, with high redshift, and the conclusion that q is negative 
remains [2]. Since then, there has been an extensive discussion on the quality 
of the data. This led to a restricted sample of 157 supernova, called the "gold 
sample" [3]. A more recent survey led to the so-called "legacy sample", of 
about 100 supernova, with high quality data [3]. Even if the precise estimation 
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of the cosmological parameters, like the raatter density, Hubble parameter, etc, 
depends quite strongly on the choice of the sample, the conclusion that the 
Universe is accelerating has remained. Hence, a large part of the community of 
cosmologist accepts the present acceleration of the Universe as a fact. 

A combination data from CMB anisotropics of the cosmic microwave back- 
ground radiation [5], large scale structure [5] and type la supernovae data [7], 
indicates an almost flat Universe, fi^ ~ 1 and a matter (zero effective pressure) 
density parameter of order flm ~ 0.3 [5]. Since an accelerated expansion can 
be driven by a repulsive effect, which can be provided by an exotic fluid with 
negative pressure, it has been concluded that the Universe is also filled by an 
exotic component, called dark energy, with density parameter flc ~ 0.7. This 
exotic fluid leads to an accelerated expansion, remaining at same time smoothly 
distributed, not appearing in the local matter clustering. 

The first natural candidate to represent dark energy is a cosmological con- 
stant, which faces, however, many well-known problems. More recently, other 
candidates have been studied in the literature: quintessence, k-essence, Chap- 
lygin gas, among many others. For a review of these proposals, see reference 
[9]. There are also claims that a phantom field (fields with a large negative 
pressure such that all energy conditions are violated) leads to the best fit of 
the observational data [TD]. A phantom field implies a singularity in a finite 
future proper time, which has been named big rip, where density and curvature 
diverge. This is of course an undesirable feature, but more detailed theoretical 
and observational analyses must be made in order to verify this scenario. Some 
authors state, based on considerations about the evaluation of the cosmological 
parameters, that there is no such phantom menace This is still an object 
of debate. 

Most of the studies made until now lay on the assumption of a simple re- 
lation between pressure and density expressed generically, in a hydrodynamical 
representation, hy p = wp°'. Quintessence, like others dark energy candidates, 
imply that w varies with the redshift, not being a constant. Chaplygin gas mod- 
els (generalized or not) |12| imply a general value for a, but typically negative, 
and w < 0. Phantom fields could be represented byQ; = l,w<— 1. Due to the 
high speculative nature of the dark energy component, many possibilities have 
been considered in the literature, both from fundamental or phenomenological 
point of views. 

In the present work we intend to exploit a more generic relation between 
pressure and density with respect to those cases normally considered. The main 
idea is to use a double component equation of state. The relation between 
pressure and density may be written as 

Pe = -kiPZ - ^2Pc • (1) 

where a, ki and k2 are constants, and the subscript c indicates that such relation 
concerns the dark energy component of the matter content of the Universe. Let 
us call the component labelled by fci as the polytropic component, and that one 
labelled by fc2 as the barotropic com,ponent. This kind of equation of state has 
been, for example, studied in a theoretical sense in reference [13j . 
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The equation of state dJ) may be also seen as a realisation of the so-called 
modified Chaplygin gas [TH [151 E] • The usual Chaplygin gas model, generalised 
or not, has been introduced in order to obtain an interpolation between a matter 
dominated era and a de Sitter phase [T2]. The modified Chaplygin gas model 
allows to obtain an interpolation between, for example, a radiative era and a 
ACDM era. In general, in this case, a negative value for a is considered. But, 
as it will be seen later, even for a positive value of a, such an interpolation is 
possible. From the fundamental point of view, the equation of state ([T]) can 
be obtained in terms of self-interacting scalar field. In references [TH [TSl [H] a 
connection with the rolling tachyon model has been established. This allows to 
consider the equation of state (IT|) as a phcnomenological realisation of a string 
specific configuration. 

In reference [17], a structure similar to ([1]) has been analysed, using obser- 
vational data, but fixing = 1, with the conclusion that the fitting of the 
supernova data are quite insensitive to the parameter a. Here, we follow an- 
other approach: we will fix a = 1/2, leaving k2 free. This has the advantage 
of leading to explicit analytical expressions for the evolution of the Universe. 
Moreover, and perhaps more important, this may lead to interesting scenarios 
where, for example, the Universe evolves asymptotically as in a de Sitter phase, 
even if the equation of state is not characteristic of the vacuum state, p = —p. 

We will test the equation of state ([T|) against type la supernovae data. We 
will span a four dimensional phase space, using as free parameters the dark 
matter density parameter QmOi the exotic fluid density parameter Qco (or al- 
ternatively, the curvature parameter flko), the Hubble parameter Hq, and the 
equation of state parameter ki or k2 . The subscript indicates that all these 
quantities are evaluated today. Using the gold sample, we will show that the 
preferred values indicate k2 ^ 6, and flco ^ 0.3, and Hq ^ 67. An open 
universe is a general prediction for this model. A phantom behaviour is largelly 
favoured. However, the dispersion is very high, and an asymptotic cosmological 
constant phase can not be discarded. 

The use of other observables, like the spectrum of the anisotropy of the cos- 
mic microwave background radiation (CMB) and the matter power spectrum, 
can in principle restrict more severely the parameter space. However, we post- 
pone this evaluation to a future study because, in both cases, a perturbative 
analysis of the model is necessary. In this case we must replace the hydrodynam- 
ical representation presented above by a fundamental description of the fluid, 
for example, in terms of self-interacting scalar fields. We note en passant that 
the hydrodynamical representation employed here may lead, at perturbative 
level, to instabilities at small scales due to a imaginary effective sound velocity, 
instabilities that can be avoided with a fundamental representation !18j . There 
are many different ways to implement this more fundamental description, which 
can lead to different results. The supernova data, on the other hand, test essen- 
tially the background, which is somehow independent of the description of the 
ffuid. 

The paper is organised as follows. In the next section, we obtain some ana- 
lytical expressions for the evolution of the Universe, and derive the luminosity 
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distance relation for the model. In section 3, we make the comparison between 
the theoretical model and the observational data. In section 4, we present our 
conclusions. 



2 The evolution of the Universe 

Let us consider the equations of motion when the exotic fluid given by the equa- 
tion of state HI) dominates the matter content of the Universe. The Friedmann's 
equation and the conservation of the energy-momentum tensor read, 



Pc 
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Pc 



SttG 







(2) 
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where k is the curvature of the spatial section. Inserting equation ([T]), with 
a = 1/2, in equation ([3]), it comes out that the exotic fluid density depends on 
the scale factor as 

fci + voa 



/32 



1/3 



(4) 



where (3 = 1 — fc2. Introducing this result in equation ([2]), it is possible to obtain 
an explicit solution for the scale factor when k = 0: 



ao< exp 



kiM t 



-Co 



(5) 



where oo and cq are integration constants. The constants obey the relations 

1/3 
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This solution can always represent an expanding Universe, with an initial sin- 
gularity. Moreover, when k2 > 1, the density goes to infinity as the scale factor 
goes to infinity, in a finite proper time, characterising a big rip. However, if 
fc2 < 1, the expansion lasts forever, and becomes asymptotically de Sitter even 
if A;2 7^ 1 (the strict cosmological constant case). Initially, the scale factor 
behaves as in the pure barotropic case with p ~ ~k2p. 

The particular case where = 1, but with free a, has been analysed in 
reference [17]. For our case, fixing a = 1/2 and k2 = I, the relation between 
density and scale factor becomes. 



-fci — 1 llna 



(7) 



If we consider the dynamics of a universe driven by the exotic fluid defined 
by equation ([T]) and pressureless matter, the equations of motion are given by, 



SttG 



(8) 



4 



Pm + = , (9) 

a 

Pc + 3-(pc+Pc) = , (10) 
a 

where pc is given by equationdl]). The conservation equations (|9|10p can be 
integrated, again for a = 1/2, leading to the relation dU and pm = Pmo/a^- In 
this case, it does not seem possible to obtain a closed expression for the scale 
factor in terms of the cosmic time t as before. However, the inclusion of the 
pressureless component is essential in order to take into account the effects of 
the baryons in the determination of the allowed range for the parameters of the 
model using the supernova data, as it will be done in the next section. 



3 Fitting type la supernovae data 

As time goes on, more and more high redshift type la supernovae are detected. 
Today, about 300 high z SN la have been reported. However, there are many 
discussions on the quality of these data. A " gold sample" , with the better SN 
la data, with a number of 157 SN, has been proposed [3]. More recently, the 
Supernova Legacy Survey (SNLS) was made public, containing around 100 SNIa 
[3]. In this work, we will use the gold sample. This will allows us to compare 
our results with previous ones using a similar method, but with different models 

From now on, we will normalise the scale factor, making it equal to one 
today: oo = 1. Hence, the relation between the scale factor a and the redshift 
z becomes \ + z = 1/ a. In order to compare the observational data with the 
theoretical values, the fundamental quantity is the luminosity distance [El [20] , 
given by 

DL = (l + z)r , (11) 

where r is the comoving radial position of the supernova. For a flat Universe, 
the comoving radial coordinate is given by 

dz' 

(12) 



H{z') 

Using equation ([H]), with the expressions for the exotic and pressureless fluid in 
terms of a, converted to relations for those components in terms of the redshift 
z, we obtain the dependence of the Hubble parameter in terms of z. Hence, the 
final expression for the luminosity distance is, for our model with ^2 ^ 1, 

Dl - ^-^^^ rin,no{i + z'f + r!fco(i + z'f 

2^ -1/2 

dz' . (13) 



Ho Jo 
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For the case ^2 = 1, the luminosity distance is given by 

1 + 
Ho 



inmQ{i + z'f + nko{i + z'Y + n,a{i-iki in(i + z') 



In the expressions above, Hq is the Hubble parameter today, which can be 
parametrized by h, such that = 100 h km/ s-Mpc. The parameter ki has been 
redefined as fci/^PcO ^i, Pco being the exotic component density today. This 
redefinition is made in order to obtain a dimensionless parameter ki. Moreover, 

rimo = 8^Gp™o/(3Ho'), ^co = 87rGpeo/(3i/o') and ^ko = -k/Hl 

The comparison with the observational data is made by computing the dis- 
tance modulus, defined as 

which is directly connected with the difference between the apparent and abso- 
lute magnitudes of the supernovae. The quality of the fitting is given by: 

where /Xq^ and /Iq^ are the observed and calculated distance moduli for the ith 
supernova, respectively, while tjQj is the error in the observational data, taking 
already into account the effect of the dispersion due to the peculiar velocity. 

In principle, the model contains five free parameters: ki, k2, Hq, O^o and 
f^mo- We will work in a four dimensional phase space: for each value of fci, we 
will vary the other four parameters. Thus, the function will depend on four 
parameters, flco, ^mOi Hq and k2. The probability distribution is then given by 

Pin,o,n„,o,Ho,k2)=Ae , (17) 

where A is a normalisation constant and a is directly related to the confidence 
region. The graphics and the parameter estimations were made using the soft- 
ware BETOCS EI]. 

The multidimensional plot of the probability distribution is not, in general, 
the best way to have an overview of the results. However, we can construct a 
two dimensional probability distribution, integrating in two of the parameters. 
In figure 1, we display the Probability Density Function, PDF, after integrating 
the distribution (jl7p in the variables Hq and k2 : a two-dimensional probability 
distribution is obtained for the variables and ficOj and with ki — and 
1.0. The plots show the confidence regions at la (68%), 2a (95%) and 3cr (99%) 
levels. This two-dimensional probability distribution reveals that the matter 
density parameters have their higher probability around flrnOj^cO ~ 0-3. An 
open model is clearly preferred. This is also evident in figure 2, where the two- 
dimensional probability distribution for flko and is shown. Such preference 
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Figure 1: The plots of the joint PDF as function of (HcO, f^mo) for the two component model, 
for fci = and 1.0. The joint PDF peak is shown by the large dot, the confidence regions 
of Icr (68,27%) by the dotted line, the 2cr (95,45%) in dashed line and the 3 a (99,73%) in 
dashed-dotted line. 
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Figure 2: The plots of the joint PDF as function of {fikO' ^mo) for the two component model, 
for fci = and 1.0. The joint PDF peak is shown by the large dot, the confidence regions 
of Icr (68,27%) by the dotted line, the 2cr (95,45%) in dashed line and the 3a (99,73%) in 
dashed-dotted line. 



for an open model contrast strongly with a similar analysis for the ACDM and 
Chaplygin gas models, for which a closed universe is clearly preferred [7]. As 
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Figure 3: The plots of the joint PDF as function of (A;2, HcO) for the two component model, 
for fci = and 1.0. The joint PDF peak is shown by the large dot, the confidence regions 
of Icr (68,27%) by the dotted line, the 2cr (95,45%) in dashed line and the Scr (99,73%) in 
dashed-dotted line. 
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Figure 4: The plots of the joint PDF as function of {k2,Ho) for the two component model, 
for fci = and 1.0. The joint PDF peak is shown by the large dot, the confidence regions 
of Icr (68,27%) by the dotted line, the 2cr (95,45%) in dashed line and the 3a (99,73%) in 
dashed-dotted line. 



the value of ki grows, a low density universe becomes more favoured, as it can 
be seen in table 1. In this table, we include also the case ki = 0.5 to show more 
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Figure 5: The plots of the PDF as function of HmO for the two component model, for fci = 
and 1.0. The joint PDF peak is shown by the large dot, the confidence regions of 1 cr (68, 27%) 
by the dotted line, the 2cr (95,45%) in dashed line and the 3 a (99,73%) in dashed-dotted 
line. 

explicitely that the parameter estimations change very slightly with ki . 

In figure 3, the two-dimensional PDF is displayed when we integrate on 
flmo a-iid Hq, remaining with k2 and f2cO- As in the preceding case, a low 
density parameter for the dark energy component is preferred. What is specially 
interesting is that the values of k2 larger than 1 and until around 40 have higher 
probabilities. This means that that a phantom scenario is clearly favoured. The 
peak of probability for k2 occurs near 5 — 6. As the value of ki increases, the 
peak probability for fc2 occurs at a lower value. 

In figure 4, the two-dimensional PDF is displayed when we integrate on Qmo 
and f^cOj obtaining a two-dimensional graphic for k2 and Hq. It is interesting 
to note, now, that regions around Hq ^ IQkm/Mpc.s are preferred. This may 
reconcile the estimations obtained using supernova with those obtained using 
CMB and matter clustering, which indicates Hq around 72km/Mpc.s [51 [^. 
The prefered value for Hq depends very little on fci. 

A more precise estimation of the parameters can be obtained by evaluating 
the one-dimensional PDF, by integrating on the three other parameters. The 
results are displayed in figures 5, 6, 7 and 8, and confirm what has been said 
based on the two-dimensional graphics. The error bars at Icr, 2a and 3a con- 
fidence levels are indicated. For the matter density parameter, its preferred 
value varies from ^Imo — 0.233 for fci = to flmo = 0.228, for fci = 1.0, with a 
very limited dispersion. The preferred value for fico decreases also from 0.347 
to 0.343. For Hq, the preferred value remains Hq ~ 67.01fcm/s • Mpc with a 
small dispersion. However, for k2, the preferred value varies from k2 — 5.75 for 
ki — to k2 ^ 4.88 for fci = 1.0. However, the dispersion is extremely large. 
Observe that fci = implies that our model reduces to a single component dark 
energy fiuid. In this case, phantom fiuids (fc2 > 1) are clearly preferred. Hence, 
an open universe, dominated by a phantom field, is a generic prediction of the 
model. Notice that the age of the universe is compatible with other astrophysical 
estimations, remaining around 13.6 Gy. 
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Fi gUre 61 The plots of the PDF as function of (llco) for the two component model, for ki = 
and 1.0. The joint PDF peak is shown by the large dot, the confidence regions of 1 cr (68, 27%) 
by the dotted line, the 2 a (95,45%) in dashed line and the 3 a (99,73%) in dashed-dotted 
line. 



PDF for DCM PDF for DCM 




Fi gure 7'. The plots of the PDF as function of {k2) for the two component model, for A:i = 
and 1.0. The joint PDF peak is shown by the large dot, the confidence regions of 1 cr (68, 27%) 
by the dotted line, the 2cr (95,45%) in dashed line and the 3cr (99,73%) in dashed-dotted 
line. 

In order to make a proper comparison, the same analysis for the ACDM, us- 
ing the same supernova sample, leads to fimo — l-Ol-o 085 ^^"^ ^0 — 65.0^174 km/ s- 
Mpc [7]. Hence, the double component fluid model exploited here predicts a low 
density universe in strong contrast with the ACDM model. The lowest value 
for the is 1-10 four the double component model, and 1.11 for the ACDM 
model. This shows that the double component model is quite competitive. 
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Figure 8: The plots of the PDF as function of (h) for the two component model, for fci = 
and 1.0. The joint PDF peak is shown by the large dot, the confidence regions of 1 (j (68, 27%) 
by the dotted line, the 2 a (95,45%) in dashed line and the 3 a (99,73%) in dashed-dotted 
line. 



ki 


0.0 


0.5 


1.0 




n 900+0. 303 
V.ZOO_Q i-^l 


n nqi+u.3Ub 


n 000+0.0. 3U4 
"•^^0-0.172 




Q 047+o.iyi 

-0.181 


^•'J^^-0.178 








Q 4.2+0.317 


(.4-, .+0.318 


Ho 


67.0li«:«^ 




67.011^2! 


k2 


t- 7r+28.V4 
^■'^-4.75 


C- qr+^8.8y 


4.88tS- 


to 




13.65li:^,^ 


13.65+^075 



Table 1: Estimated values for fl„io, flcO, flko, Hq, ^2 and to for three different 
values of fci, at 2a level. 

4 Conclusion 

In this work, we have explored the possibility that the dark energy has an 
equation of state given by ([T]). This proposal has already been explored in ref- 
erence [T7] , but restricting one of the components to behave like a cosmological 
constant: a two-component fluid, which is a "variation" around a cosmological 
constant, has been analysed. Here, we alleviate this restriction, but introduc- 
ing another one: the linear component can have any barotropic index, but the 
second component must vary as the square root of the density. This allows us 
to obtain an analytical expression for the evolution of the Universe, at least for 
a flat spatial section. This analytic expression reveals that it is possible to have 
an asymptotically de Sitter phase, for p < and ^2 < 1, even if only ^2 = 1 
represents the cosmological constant. This is an intriguing aspect of the model. 
For fc2 > 1 there is always a big rip. 

The restriction in the polytropic factor a seems not to be so relevant in 
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view of the results of reference [T7] : if the second component obeys a polytropic 
power law, there is a strong degeneracy on the polytropic factor, and almost 
any value of the power is allowed. 

Our results indicate an open universe with a matter density parameter 
around ^Imo ^ 0.3, with a similar estimation for the dark energy component. 
The ACDM model favours a closed universe. On the other hand, the predicted 
value for the Hubble parameter is more consistent with other observational tests, 
like CMB, that is, Hq ~ 67 km/ s ■ Mpc [53]. For the barotropic index in the 
two-component fluid the results indicate that /c2 > 1 is highly favoured. The 
dispersion is very high, but tends still to favour a phantom scenario. As the 
component index fci increases, the preferred value of fc2 decreases slightly. 

It must be remarked that the model described here exhibits a slightly 
smaller than for the ACDM model. This shows that the model is quite com- 
petitive. In our opinion, even if a phantom scenario is clearly preferred, the 
fact of predicting a low density universe is interesting in its own, mainly when 
compared with the analysis of clustering of matter in the universe. 
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